Finite temperature Casimir effect for massless Majorana fermions in a magnetic field 



I. INTRODUCTION 
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The zeta function regularization technique is used to study the finite temperature Casimir effect 
for a massless Majorana fermion field confined between parallel plates and satisfying bag boundary 
conditions. A magnetic field perpendicular to the plates is included. An expression for the zeta 
function is obtained, which is exact to all orders in the magnetic field strength, temperature and 
plate distance. The zeta function is used to calculate the Helmholtz free energy of the Majorana 
field and the pressure on the plates, in the case of weak magnetic field and strong magnetic field. 
In both cases, simple analytic expressions are obtained for the free energy and pressure which are 
very accurate and valid for all values of the temperature and plate distance. 
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The Casimir effect is a purely quantum phenomena where an attractive or repulsive force occurs between electrically 
neutral conducting plates in vacuum. The effect was first discovered by Casimir who calculated the attractive 
electromagnetic force between two parallel conducting plates, and can be regarded as a quantitative proof of the 
quantum fluctuations of the electromagnetic field. The repulsive Casimir effect was discovered by Boyer some time 
Q-T later Q, when he showed that if the electromagnetic field is confined inside a perfectly conducting sphere, the wall 
of the sphere is subject to a repulsive force. The first experimental evidence of the Casimir force was obtained more 
than fifty years ago by Sparnaay |!3j and, since then, many greatly improved experimental observations of this effect 
have been reported. For a comprehensive review of these experiments see the review article by Bordag et al. [4|. 
CN ' Since Casimir forces have many applications, from nanotubes and nanotechnology d, @ , to branes and compactified 
extra dimensions @-Hl] and to string theory 0, [Hj] , a great deal of effort has gone into studying the Casimir effect 
and its generalization to quantum fields other than the electromagnetic: fermions, bosons and other scalar fields have 
| all been investigated extensively Q. 

. It is well known that the Casimir effect is very sensitive to the boundary conditions for all types of quantum fields, 
and the most used ones are Dirichlet and Neuman boundary conditions on the plates. These boundary conditions 
however cannot be used in the case of fermion fields or fields with spin in general (26| and therefore, for fermions, the 
bag boundary conditions are used, which originally were introduced to find a solution to confinement [13] ■ I n this 
work I will use the bag boundary conditions for a Majorana fermion field confined between two parallel plates. 
Majorana fermions appear in many different areas of physics: in some models of extra dimensions a Majorana 
.£h ! fermion, the Kaluza-Klein neutrino, is a leading WIMP candidate [2S| . In the field of superconductivity, a Majorana 
bound state is theoretically predicted in rotating superfluid 3 He - A between parallel plates and in the presence of 
a magnetic field [2t| . While the Casimir effect for Majorana fermion fields has been studied in vacuum (3^, [HI and 
recently at finite temperature [32j , a study that considers magnetic effects in vacuum and thermal and magnetic effects 
combined has not appeared in the literature. This paper will therefore investigate the Casimir effect for Majorana 
fermion fields at finite temperature and in the presence of a magnetic field. 

Casimir effect calculations must address carefully the issue of the vacuum energy regularization. Many regularization 
techniques are available nowadays, and two of them have been applied frequently and successfully to the Casimir effect: 
the cutoff method often used in various piston configurations (3lL |33| and the zeta function technique [34- 36] . My 
choice for this paper is the zeta function technique, which is very powerful and is also used in the computation of 
effective actions [37l . |38| . I will apply this regularization method to calculate the free energy and pressure on the 
plates for a Majorana fermion field with bag boundary conditions confined between two parallel plates, at a distance 
a from each other. I assume this system to be in thermal equilibrium with a heat reservoir at finite temperature T, 
and will use the imaginary time formalism of finite temperature field theory, which is suitable for a system in thermal 
equilibrium. A magnetic field B is present in the region between the plates and is perpendicular to the plates. 
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In Section HH I calculate the zeta function for this system, exact to all orders in eB, T and a, where e is the Majorana 
field charge. In Section IIIIl I use the zeta function obtained in the previous section, to calculate the Helmholtz free 
energy of the Majorana fermion field and the pressure on the plates and obtain simple analytic expressions for these 
quantities in the case of weak magnetic field (eB -C T 2 ,a~ 2 ), and of strong magnetic field (eB ^> T 2 ,a~ 2 ). A 
discussion of my results is presented in Section HVl 



II. ZETA FUNCTION EVALUATION 



Using the imaginary time formalism of finite temperature field theory I write the partition function Z for a fermionic 
system in thermal equilibrium at finite temperature T 



N I D4> Dip exp [ f dr j d 3 xC\ (1) 

J Antipcriodic \J0 J I 



where C is the Lagrangian density for the fermionic system, N is a constant and 'antipcriodic' means that this 
functional integral is evaluated over field configurations satisfying 

ip(x,y,z,0) = -i/)(x,y,z,P), (2) 

where j3 = 1/T is the antipcriodic length in the Euclidean time axis. In addition to the boundary conditions given by 
@, I impose bag boundary conditions [271 ] for Majorana fermions between two plates. In 3-dimensional space with 
two very large parallel plates perpendicular to the z-axis and located at z = and z — a, the MIT bag boundary 
conditions constrain the Majorana fermion fields to the region between the plates. These boundary conditions are 
expressed in covariant form as 

n^^rp = (3) 

where = (0, ft) and ft is the vector normal to the surface of the plates and directed towards the interior of the slab 
region. The above boundary conditions prevent the flow of fermion current out of the slab region. In the slab region 
there is also a uniform magnetic field pointing in the z direction, B = (0,0, B). The Majorana fields have charge e 
and will interact with the magnetic field. 

The Helmholtz free energy F and the partition function Z are related by 

F = -r 1 logZ, (4) 

therefore I evaluate the functional integral ((T|) and, after the straightforward gamma algebra, find 

2 log Z = log det (-D^l-Fa) + log det \T a ) (5) 

where the symbol T a indicates the set of functions which satisfy boundary conditions ([2]) and (|3]), and the operators 
DjQ are defined as: 

D^=d* + d 2 z -(p-eA)i T eB, (6) 
where A is the electromagnetic vector potential, the subscript E indicates Euclidean time, and I used the notation 

P± = (Px,Py,0). 

The zeta function technique allows me to use the eigenvalues of — to evaluate \ogZ. The bag boundary 
conditions (J3|) are satisfied only if the allowed values for the momentum in the z-direction are 

where I £ {0,1,2,3,...}, and therefore the eigenvalues of — <9 2 — 9 2 whose eigenfunctions satisfy and ([3]) are: 

7T 2 / 1\ 2 4^ 2 / 
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where I e {0,1,2,3,...} and m € {0, ±1, ±2, ±3, ...}. The spectrum of the operator (p— eA)\ is well known from 
one-particle quantum mechanics, and its eigenvalues are 



2eB n 



1 



(9) 



with n G {0,1,2,3,...}. Using the eigenvalues (JH]) and (j9]), I construct the zeta functions C,(s,— D^j and 
C f s, ~£>e ') ' w hich are given by 



oo oo 



/— m=— oo x 7 n— 

where L 2 is the area of the plates. I add them to obtain 

c( s ) = c( s ,-4 +) )+ c ( s '^ } )' 

and find immediately, using the zeta function technique 

21ogZ = -C'(0). 

Using the following identity 

1 



7T /, 1\ 2 4tt 2 



/? 2 



m+~) + 2eB\n+-)±eB 



r(s) Jo 



I rewrite C(s) as 



C(«) = 



where I also used 



L 2 1 
27rr(s) 



EE/ dtt s ~ 2 eBt coth eBtexp 

] — n to — — rsr> *^ 



/— m— — oo 



7.- /, 1\ 2 I"" /' 1 V ' 



B 2 



e 2nz + e - 2 ("+J> 



= coth ; 



After applying the Poisson resummation formula 39] to the two sums, I obtain 

V 13 1 dtt s ~ 3 eBt coth eBt '--tVi-i,', 1 \ I i ^ i ><", - " 



C(«) 



2^ 2 r(s) 7 



fi+Ec-D'e-'^ (Ve 

V 1=1 / V m=i 



(-l) m e- 



(10) 

(11) 
(12) 

(13) 
(14) 
(15) 

(16) 



where V — L 2 a is the volume of the slab. It is evident from fp~6|) that the Poisson resummation formula allows me to 
naturally divide £ into four parts 



COO = Che,o(s) + Che,p(s) + Cc,o(s) + (c,p(s), 



(17) 



two of which are vacuum parts that do not depend on (3 except for the overall multiplicative factor, while the other 
two are finite temperature corrections that vanish when j3 — > oo. The four parts of £(s) are: 



Che,o(s) 



L 2 aj3 



dtt s ~ 3 eBt cotheSt, 



L 2 aP 
4tt 2 r(s) 



/ dtt s - 3 eBt coth eBt V(-l) 

J ° Vm=l 



(18) 



(19) 



Cc,o(s) 



r2 a R r°° I °° 

— tt—tt dtt s ~ 3 eBt coth eBt V 
4tt 2 T( S ) y ^ 



(_l)i e -« a «V* , 



(20) 
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C ^ (S) = dtt^eBtcotheBt [^£{-1)^ >> j (^(-1)^^/^ , (21) 

and each one will contribute differently to the free energy. I will show in the next section that the contribution of 
(he,o to the free energy is the opposite of the unrenormalized Heisenberg-Euler effective lagrangian [4l| multiplied by 
the volume of the slab V, while the contribution of Che,0 is the opposite of the finite temperature correction to the 
Heisenberg-Euler effective lagrangian, also multiplied by V. The two remaining pieces, Ceo and Cc,/3i will be shown 
to be the main contributors to the vacuum Casimir energy and its finite temperature correction. 

III. FREE ENERGY AND PRESSURE 

I combine ^ and (fT2"|) to obtain the free energy F in terms of the zeta function 

F =jpC(0). (22) 

The derivatives of the four parts of the zeta function are obtained easily by taking advantage of the useful fact 
that, for a well behaved G(s), the derivative of G(s)/T(s) at s = is simply G(0). I find 

Ch£o(°) = f4/ dtt- 3 eBtcotheBt, (23) 



Stt 2 Jo 

and its contribution to the free energy is 



Fhe,o = t^I dtt~ a eBtcotheBt, (24) 







where we recognize immediately that Fhe,o — —V£he,0i where £he,o is the unrenormalized Heisenberg-Euler 
effective lagrangian for massless Majorana fcrmions. Similarly I find 

Che A ) = r dtr 3 eBtcotheBt ( £ (_i)™ e -^ 2 ™ 2 /« j ; (25) 

whose contribution to the free energy is 

V f°° ( i 

Fhe,P = J~2 dtt~ 3 eBt coth eBt ^ (_ x yn e ~ t> 2 m 2 / it ( ( 26 ) 

\m— 1 / 

and Fhe,p = — V£he,/3, where Che,/3 is the finite temperature correction to the Heisenberg-Euler effective lagrangian 
[42l | for massless Majorana fermions. Notice that Fhe,o and Fhe,p do not depend on the plates distance a except for 
a multiplicative factor of a contained in V. The remaining two parts of the zeta function, C,c,o and Cc,p, are the ones 
mostly contributing to the Casimir force, and their derivatives are given by 



C^ o (0) = ^/ dtt- 3 eBtcotheBt[J2(-l) l e- a2l2 / t ) (27) 



J dtt~ 3 eBt coth eBt ( ^(-l)'e~ a2 ' 2 /* J 



and 



2tt 2 

Their contributions to the free energy are 



CcA Q ) = ^l°°dtt-*eBtcotheBt ^(-lj'e-" 1 *^ (f, (-l)- e -^/« ) . (28) 



2/3 

and 



^,0 = ^# (29) 



*b,/J = (30) 



5 



respectively. It is not possible to evaluate (128)) in closed form for arbitrary values of B, a and (3, but it is possible 
to find simple expressions for (|27|l and (|2"5)) when the magnetic field is weak, eB <C a~ 2 , f3~ 2 , and when B is strong, 
eB > a- 2 ,(3~ 2 . 

In the case of weak magnetic field I can set 



eBtcotheBt w 1 + Tr(eBt) 2 
o 

and, after substituting (f5T|) into (|2"T)) . integrate by changing variable from i to -, to find 



(31) 



(32) 



where Cr(4) = 7r 4 /90 is the Riemann zeta function of number theory, T(2) = 1 is the Euler gamma function and 



I find 



f(s,a) = a 2s (2 



df(s = 0,a) 
ds 



2s/ 2s+l 



l)Cfl(-2«) 



£(-£) 

r( s ) ■ 



7£ + 2C^(0)+ln(4a), 



(33) 



(34) 



where je — 0.5772 is the EulerMascheroni constant and Cfi(0) = —0.9189 is the derivative of the Riemann zeta 
function. The interesting numerical fact 



allows me to write 



7£ + 2^(0)+ln4 
df(s = 0,a 



ds 



= In a H — , 



and, after inserting ([36]) into J32J) and then ((32} into (|29j>. I find 

,'7^ tt 2 V , (eB) 2 



8 J 720 a 4 24tt 2 



V In a 



(35) 



(36) 



(37) 



which, for B = 0, is the vacuum Casimir energy for massless Majorana fermions with bag boundary conditions and is 
in agreement with (43[, and with |3lj . where the same result is obtained by a different method, using a piston model 
to derive the vacuum Casimir energy for Majorana fermions. Notice that the second term in (|3T[) is the correction to 
the Casimir energy in vacuum due to the weak magnetic field. 

In the case of weak magnetic field I insert ([31]) into (|28[) and integrate as described previously, to obtain 



h(a,/3)T(2) 



(eB) 2 dg(s = 0,a,P) 



ds 



where the functions g(s,a,j3) and h(a,/3) are given by 
g(s,a,/3) 



(1 + A S+1 )E 2 ( -a; a 2 , £ ) - 2E 2 (-a; a 2 , /3 2 ) - 2E 2 ( -s; 4a 2 , tL. 



P 2 \] T(-s) 
4 I T{s) ' 



(38) 



(39) 



h{a,P) = \e 2 (2; a 2 , 2E 2 (2; a 2 ,/? 2 ) - 2E 2 (2; 4a 2 , ^) 



(40) 



and I have expressed the double sum in terms of Epstein functions [36|, |4CJ, |44j which for any positive integer N are 
defined by 



00 00 



E N (s; 01,02, a N ) = ^ .... ^ 



1 1 (ainf + a 2 n2 

ni — 1 U2 — 1 TIN — 1 



a N n 2 N ) s 



(41) 



6 



The Epstein functions can be analytically continued to meromorphic functions in the complex plane [36|, |40j , and for 
N = 2 their analytic continuation is given by: 



E 2 (s; ai,a 2 ) 



i nrr( s -±) 



2 2 V a 2 T(s) 



1 



2n s 



E 1 ( S --;a 1 )+ 



J^ySjQn n,m— 1 



2irmn 



(42) 



where K v {z) are modified Bessel functions. This analytic continuation of the Epstein functions allows me to find 



h(a,/3) 



7tt 4 1 



2n 2 



90 £ 4 a §pi ^ 

" n,m— 1 

and 

dg(s = 0,a,P) SPY 



^ (to) 



if a 



9s 



\ a ) ^ (to) 



V2 



V2 



TTTYlTl/3 



Ki 



2a 



7rmnf3 
2a 



2Ki 



2Ki 



2-Kmnfi 



I 2-Kmnj3 \ 
\ a j 



(43) 



In 



where I used (|35[) . With the help of the following 



(n + fc)! 



2z k\{n - k)\{2z) k 
for the modified Bessel functions of half-integral order, I find 



(44) 
(45) 



^^ = -90-^ + ^3 L ^3 



2 (1 + 



irnmL 



nnm(3\ »i-j< / 1 Trnmfi 



and 



9g(s = 0,a,/3) 
9s 



= 2 E 1 

fa * TO 



2 6 



— e a 



1 


KM 


~ 2 





(46) 
(47) 



From (|46p and (|4"7)l I can easily obtain the low temperature expression of /i and For aT < 1, I find 



/i(a,/?) 



7tt 4 1 



2tt 2 ( 2a 



and 



9g(s = 0,a,/3) 
9s 



-2e^ - i 
2 



7T/3 



In 



_i2 



(48) 



(49) 



I insert and (|4"9"|) into (f3"5)l . then (|3"5)l into (J3UJ) and find immediately the finite temperature correction to the 
Casimir energy in the case of weak magnetic field and low temperature. For eB < T 2 < a~ 2 , I obtain 



_ 7n 2 V 



V 



360 /3 4 2a 2 f3 2 



1 2a \ (e£) 2 

1 -\- — — ) 6 ~ - _ , „ V 



n(3 



24tt 2 



4. — In ( I 



(50) 



In order to calculate the Hclmholtz free energy, I need to add the contribution of the Heisenberg-Euler effective 
lagrangian (|24l) and of its finite temperature correction (f26|) to the Casimir energy and its thermal correction. In the 
case of weak magnetic field, eB <C T 2 , these contributions are [3^, H^, 



HE '° + HE ' 360 /3 4 + 48tt 2 



In 



(3 2 eB 



(51) 



and, once I add (|5ip to the vacuum Casimir energy for weak magnetic field (|37p and to its finite temperature correction 
(f50|) . I obtain the Helmholtz free energy F in the limit eB -C T 2 <C a -2 



r 



8 7 720 a 4 2a 2 /3 2 



2a \ {eB 



7T/3 7 



24tt 2 



l£ 1 

2 



4 e -* _ -lnfeBa 2 ) 



(52) 
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Notice the cancellation of the Stefan-Boltzmann term and of the temperature dependence from the logarithmic term. 
Due to these cancellations, I find that the leading order correction to the Casimir vacuum energy is the temperature 
independent and magnetic field dependent logarithmic term which, even in the case of weak magnetic field, is larger 
than the exponentially suppressed part, for Ta <C 1. The pressure P on the plates is given by 



P 



1 dF 



(53) 



and therefore, in the limit eB < T 2 « n 2 , I find 



8 J 240a 4 4a 3 8 



7T (eB) 
e ' 



12tt 



— e ^ - -hx(eBar) 



(54) 



Another interesting limit, for the weak field case, is the high temperature case of eB <C a~ 2 <C T 2 . Exploiting the 
fact that (|28p is invariant if I exchange a and /3/2, I write h and g in a different form, which is better suited for the 
high temperature expansion, aT ^> 1: 



h(a,8) 



7\ n^J_ 

8J 180 a 4 ' 2a 3 8 



7T x 1 

T77 ° 



In 



4irnma\ 



Aimma\ _2* 



& ) 



1 A%nma\ _& jnls _ 
- H e -9 

2 8 



(55) 



and 



dg(s = 0,a,p) 
ds 



-2 v -n 



n,m— 1 



m V 2 



- I In a i - 



(56) 



from which I find immediately the high temperature expansions of h and — 



ds 



h(a,P) 



7tt 4 1 



and 



1,440 a 4 a 3 8 
dg(s = 0,a,8) 



ds 



-2e" 



27ra 



- - | In a + - 



(57) 



(58) 



I use the last two equations to obtain the high temperature limit of F ( 



7 \ 7T V IV 



8 j 720 a 4 2 a 2 ^ 2 



M (eB) 



27ra 7 



24tt 2 



, _2™ 1 

V 4e 3 + In a + - 



(59) 



Adding (|59"|) to (I3T1) and (IBTj) , I obtain the free energy in the weak field and high temperature limit eB <C a 2 «T 2 



7\ 7T 2 V IV 



F = - ^ ^^-^ 



8/ 45 /3 4 2 a 2 /? 2 



1 + JLV»-I 



eB) 



27ra y 



24tT< 



■V 



4e" 



2 V 4 



(60) 



Notice that the leading term here is the Stefan-Boltzmann term. The pressure in the weak field and high temperature 
limit is obtained immediately from (|60[) 



7T _2jra (eB) 2 

Sj 45 B 4 a/3 3 " 6 13 ~ 24tt 2 



8a _2 2L a 1 , 



B 2 eB 



1 



(61) 



where the dominant piece is clearly the outward pressure originating from the Stefan-Boltzmann term of the free 
energy. 

In the case of strong magnetic field 



eBtcotheBt w eBt 



(62) 
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and, after substituting ([62]) into ([27]) and integrating, I find 

Cc,o(0) = -^r(i)c*(2). 

From (|63|) I obtain immediately the vacuum piece of the Casimir energy for strong magnetic field 

_ VeB 
Fc -°-~96^- 

In order to find the thermal correction to Fc,o, I substitute (|62j) into ([28]). integrate and obtain 



E 2 1; a z , - £ 2 1; < ^ - £ 2 1; 4a^, ^ r 



/3 5 



(63) 



(64) 



(65) 



which, using the analytic continuation (|42|) of the Epstein functions and the expression of the modified Bessel functions 
(fl5"j) . can be written as 



C^(0) = VeB 



I 2 

12/3 + TTfl 



E 1 

^ ' 777. 



From (p36f I obtain Fc,p , the thermal correction to the Casimir energy in the case of strong magnetic field 



VeB VeB 

F C ,0 = 7777^ + 



2A0 2 tt&B z — ' TYl 

n.m— 1 



oo 1 

E 1 



r 



(66) 



(67) 



(68) 



I now use this equation to obtain the low temperature limit of Fc,p and, for eB 3> a 2 3> T 2 , find 

_ VeB VeB _^ 

In the case of strong magnetic field, eB ^> T 2 , the contribution of the Heisenberg-Euler effective lagrangian and of 
its finite temperature correction to the free energy is [H, HiJ |46[ 

F HE , + F HE , p = ~ R 4- V -^S- ln (feB) ■ (69) 



I add 



24 2 48tt 2 

and (f6"4")l and obtain the free energy F and the pressure P, in the limit eB ^> a -2 3> T 2 



F = 



ln ^ eB )-96^ 



V eB VeB 



7TO/3 



P = 



V{eBf 
48tt 2 

( eB ) 2 , /o2 n \ 1 eB eB 

- V In (0 eB) 77- -\ 

48tt 2 V ^ ; 96 a 2 2a 2 



(70) 



(71) 



In both of these equations, the dominant term is quadratic in the magnetic field and comes from the Heisenberg-Euler 
effective lagrangian. 

Finally I want to find F and P in the case of strong magnetic field and high temperature, eB ^> T 2 ^> a~ 2 . I again 
exploit the symmetry of (|28|) for exchange of a and 0/2 to write (|67|) in an equivalent form, which is better suited for 
expansion in the case of T 2 ^> a~ 2 



= 



VeB VeB 



96a 2 



oo 1 

E 1 



and its high temperature limit is 



iraB z — ' m 

n.m—l 



VeB VeB _ 
F c ,p = ttt-^t -e 



1 

2 6 



96a 2 



na0 



(72) 



(73) 



By adding (|73p to (fB"8")) and I obtain the free energy and the pressure in the case of strong magnetic field and 
high temperature 



F 



V eB VeB 



24 2 ira0 



P 



VjeBf 

487T 2 

i eB ? , to% n \ 1 eB 2eB - 



-e—r. 



(74) 
(75) 
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IV. DISCUSSION AND CONCLUSIONS 



In this paper I used the zeta function regularization technique to study the finite temperature Casimir effect of a 
massless Majorana fermion field confined between parallel plates and in the presence of a magnetic field perpendicular 
to the plates. I have obtained an expression for the zeta function (|16|) which is exact to all orders in the magnetic 
field strength B, and from it I have derived expressions for the Helmholtz free energy and for the pressure on the 
plates in the case of weak magnetic field (eB <C a~ 2 , T 2 ) and in the case of strong magnetic field (eB ^> a -2 , T 2 ). 

In the case of a weak magnetic field, I found the temperature correction to the Casimir energy to be 



■ C,p 



V 
4tt 2 



h{a,P) 



{eB) 2 dg{s = O,a,0) 



ds 



(76) 



and obtained two equivalent expressions for h, (|46|) and (|55j) . which are exact to all orders in a and /3, and two 
equivalent expressions for ||, ([47|) and (|56|) . which are also exact to all orders in a and (3. These expressions involve 
two infinite double sums, and I have been able to evaluate them numerically with very high precision. I find that, for 
< aT < i the simple low temperature expression of h that I obtained in Eq. (1481) is within less than one percent 
of the exact value of h, as given by the infinite double sum of Eq. (|46)) . Similarly I find that, for h < aT < oo, the 
high temperature expression of h that I write in Eq. (|57l) is within less than one percent of the exact value of h. A 
similar numerical evaluation of the exact expression of ^ leads me to discover that the low temperature expression 
of (|49|) . is within less than four percent of its exact value in the range < aT < |, while the high temperature 
expression of (|58p . is within less than four percent of its exact value in the range g < aT < oo. I summarize these 
findings by writing the free energy F for weak magnetic field as 



F = 



(1\ _zrl v — 

\8/ 720 a? 
\8j 45 4 



V 

2a 2 /3 2 

V 
2a 2 p 2 
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2a 

7T/3 
P 

2na 



-Z£. (eB) 2 
J 24tt 2 



V 



V 



Ae 



\ In (eBa 2 ) 



In 



for < aT < \ ; 
for \ < aT < oo 



(77) 



When evaluated at aT — i , both expressions of F yield the same value. Equation ([77|) is a simple analytic expression 
of F for weak magnetic field, valid for all values of the temperature T and the plate distance a, and with a discrepancy 
of no more than a few percent from the exact value of F. A similarly accurate expression of the pressure P, valid 
for weak magnetic field and all values of a and T, is obtained immediately from (|77jl. since P = —j2^;- Notice 
that, if we set B = in (|77|) . we obtain the free energy for the finite temperature Casimir effect of Majorana fermion 
fields and, opposite to what is claimed in 32[, the free energy is always negative and does not become positive for 



aT > 0.37. Even if we do not include the contribution of the thermal part of the effective Lagrangian and consider 
only Fc,o + Fc,/3 at B = 0, we find 



Fc,o + Fc.b — 



(1) 



jrf_ V 

720 7? 



+ (i)f 



v 

45 ~W 



2a 2 /3 2 



1 



2n 

nil 



for < aT < \ 



v 

' 2a 2 p 2 



P 
2na 



for 



i < aT < oo 



(78) 



which is negative for < aT < oo. 

In the case of a strong magnetic field, two equivalent expressions of the exact temperature correction to the Casimir 
energy are presented in Eqs. (|67[) and (1721) . both involving two infinite double sums. After a highly precise numerical 
evaluation of these two expressions, I find that the low temperature limit of Fc,p presented in (1681) . is accurate within 



less than two percent for < aT < |, while the high temperature limit of Fc t p shown in (fT3")l is also accurate within 
less than two percent for | < aT < oo. These findings allow me to write the free energy for strong magnetic field as 



F = 



B)' 

4Stt 2 



In ((3 2 eB) - ^ _ VeB 

\" ) 96 a 2 Trap 

v[eB T \n(3 2 cB) v eB VeB ~- 



P 2 



ir a ft 



for < aT < \ ; 
for i < aT < oo 



(79) 



a simple analytic expression of F, valid for all values of T and a, and with a discrepancy of no more than one or two 
percent from the exact value of F. The pressure in the case of strong magnetic field is obtained immediately from 
(fT9"|) for all values of a and T. 
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